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Abstract 

With the aim of a further investigation of the nonperturbative Hamiltonian approach 
in gauge field theories, the mass spectrum of QED-2 is calculated numerically by using 
the corrected Hamiltonian that was constructed previously for this theory on the light 
front. The calculations are performed for a wide range of the ratio of the fermion mass 
to the fermion charge at all values of the parameter 9 related to the vacuum angle 9. The 
results obtained in this way are compared with the results of known numerical calculations 
on a lattice in Lorentz coordinates. A method is proposed for extrapolating the values 
obtained within the infrared-regularized theory to the limit where the regularization is 
removed. The resulting spectrum agrees well with the known results in the case of 9 = 0; 
in the case of 6 = ir, there is agreement at small values of the fermion mass (below the 
phase-transition point). 
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1. Introduction 



The Hamiltonian approach to quantum-field theory in light-front coordinates x ± = (x° ± 
x 3 )/\/2, x 1 - = (a; 1 ,^ 2 ), where x + plays the role of time [1], is one of the nonperturbative meth- 
ods for solving the problem of strong interaction [2, 3]. Within this approach, the quantization 
is performed in the x + = plane, the generator P + of a shift along the x + axis playing the 
role of the Hamiltonian. The generator P_ of a shift along the x~ axis does not displace the 
quantization surface; therefore, it is a kinematical generator (according to Dirac's terminology) 
in contract to the dynamical generator P + . As a result, the momentum operator P_ appears 
to be quadratic in fields and does not depend on interaction. At the same time, the operator 
P_ is nonnegative and has zero eigenvalue only on the physical vacuum. This results in that 
the field Fourier modes corresponding to positive and negative values of p_ play the role of, 
respectively, creation and annihilation operators over the physical vacuum and can be used 
to construct Fock space. Thus, we see that, in light-front coordinates, the physical vacuum 
formally coincides with the mathematical vacuum. 

The spectrum of bound states in the theory can be sought by solving the Schroedinger 
equation 

P + |*)=p + |*> (1) 

in the subspace specified by fixed p- and p± and by employing the expression m 2 = 2p+p~ —p 2 ± 
for the mass. This search for bound states can be performed beyond perturbation theory - for 
example, with the aid of so-called method of discrete light-cone quantization [2, 4]. 

However, the light-front Hamiltonian formalism involves a specific divergence at p_ = 
[2, 3], and it must be regularized. The introduction of a cutoff > e > 0, which violates 
Lorentz and gauge invariance, is one of the methods for its regularization. A cutoff \x~ \ < L that 
involves imposing (anti)periodic boundary conditions in x~ (discrete light-cone quantization 
method, which respects gauge invariance) is yet another possible regularization. In this case, 
the lightlike momentum p_ becomes discrete (p„ = p n = nn/L, where n is an integer), the 
field zero mode corresponding to n = being separated explicitly. In principle, the canonical 
formalism makes it possible to express this zero mode in terms of other modes by solving the 
constraint equation, but this is difficult as a rule [5, 6]. 

The regularization of the above divergence usually renders a theory in light-front coordinates 
nonequivalent to its conventional formulation in Lorentz coordinates [7-9] . This can be revealed 
even in the lowest orders of perturbation theory [10]. As a result, there arises the problem of 
correcting the canonical light-front Hamiltonian (which is the result of a "naive" canonical 
quantization in light-front coordinates) - that is, the problem of seeking counterterms to it 
that compensate for the above distinctions between the Hamiltonians. If this problem can be 
solved for a specific theory in all orders of perturbation theory, the resulting corrected light-front 
Hamiltonian can then be used to perform nonperturbative calculations. 

The aforementioned formal coincidence of the physical and the mathematical vacuum be- 
comes rigorous only after the introduction of regularization, upon which the vicinity of the 
point p- = is eliminated - for example, after the introduction of the cutoff > e > 0. If 
the corrected light-front Hamiltonian can be constructed for a theory regularized in this way, 
then vacuum effects inherent in the original theory in Lorentz coordinates must be taken into 
account with the aid of additional terms of this Hamiltonian. 



2 



The problem of constructing the corrected canonical light-front Hamiltonian was success- 
fully solved both for nongauge field theories of the Yukawa model type [11] and for QCD in 
the gauge A_ = [12]. In the last case, however, the corrected light-front Hamiltonian was 
constructed only for specific ultraviolet and infrared regularizations violating gauge invariance. 
As a result, it appears that the corrected Hamiltonian involves a large number of indeterminate 
coefficients; only for some, a priori unknown, dependence of these coefficients on the regulariza- 
tion parameter does it reproduce, in the limit where the regularization is removed, the results 
of the Lorentz-covariant theory in all orders of perturbation theory. The practical calculations 
with the resulting Hamiltonian are very cumbersome because of the presence of unknown coef- 
ficients and because of a complicated structure of regularization (the regularized Hamiltonian 
involves a large number of additional fields). 

In view of these circumstances, it is desirable to seek alternative methods for constructing 
the correct light-front Hamiltonian for gauge theories. In this connection, it is of interest 
to study the simplest models that admit a nonperturbative approach - in particular, those 
where one can study the behavior of infinite series of perturbation theory in all orders. Two- 
dimensional QED (QED-2) featuring a nonzero fermion mass (it is also known as the massive 
Schwinger model) is one of such models. In recent years, this two-dimensional model has 
attracted attention as an object of application of new methods for studying QCD, since it 
possesses many properties similar to those of QCD: confinement, chiral-symmetry breaking, 
and a topological 9 vacuum (see [13] and references therein, as well as [14-16]). Information 
obtained in analyzing QED-2 can also be used in developing new methods that take into 
account nonperturbative vacuum effects and which are appropriate for constructing the light- 
front Hamiltonian for four- dimensional gauge theories. It should be noted that attempts at 
extracting information about four-dimensional gauge theories on the light front from an analysis 
of QED-2 were undertaken earlier [17]. 

For QED-2, there exists the possibility of going over to an equivalent scalar theory [18] 
(belonging to the type of the sine-Gordon model). This can be done by means of the bosoniza- 
tion procedure - that is, by going over from the fermion variables to boson ones [9, 19]. Upon 
this transition, the mass term of the fermion field in the QED-2 Hamiltonian becomes the 
interaction term for a scalar field, while the fermion mass M becomes the interaction constant 
in the boson theory. In the boson theory, the fact that the quantum vacuum in QED-2 has 
a nontrivial character associated with instantons (9 vacuum) [18, 15] is taken explicitly into 
account with the aid of the parameter 9 in the interaction term. At M — 0, QED-2 reduces to 
the Schwinger model, while the equivalent boson theory appears to be free. 

Perturbation theory for a boson theory (perturbation theory in the fermion mass) is usually 
referred to as chiral perturbation theory. For this kind of perturbation theory, ultraviolet 
finiteness was proven in [20, 21]. By analyzing perturbation theory in all orders of M, one 
can construct a corrected light-front Hamiltonian in terms of bosons and, after this, return 
to the fermion variables [21, 22]. It should be noted that boson perturbation theory differs 
radically from perturbation theory in the coupling constant of the original theory involving 
fermions (in QED-2, the latter perturbation theory does not exist at all because of infrared 
divergences, and this was the reason for introducing bosonization). Therefore, the resulting 
light-front Hamiltonian can take into account nonperturbative (in the conventional coupling 
constant) effects. But at the same time, it can fail to describe effects that are nonperturbative 
in the fermion mass - for example, phase transitions. 
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It is well known that, at least at the vacuum-angle value of 9 = it, there is a phase transition 
in QED-2 at some value of the fermion mass M (see, for example, [13]). It should be expected 
that, in the presence of a phase transition, which is accompanied by the appearance of nonzero 
vacuum expectation values of some operators, the correct light-front Hamiltonian must have 
different form for different phases, since the light-front vacuum itself is always trivial. Therefore, 
the results of calculations performed with a specific Hamiltonian must be valid only within one 
phase. The calculations performed in the present study corroborate these considerations and 
make it possible to determine the presence and an approximated position of the phase-transition 
point. A similar phenomenon was discovered previously in the simple two-dimensional Ay? 4 
scalar-field model (see [9, 23], where an approximate method was proposed for going over to 
the Hamiltonian describing a different phase). 

In this study, we perform a numerical nonperturbative calculation of the mass spectrum 
of the corrected light-front QED-2 Hamiltonian constructed in [21, 22]. The results obtained 
in this way are compared with the results of numerical calculations on a lattice in Lorentz 
coordinates from [14, 13]. 



2. Corrected light-front Hamiltonian for QED-2 

The QED-2 Lagrangian density in Lorentz coordinates has the form 

L = --F^ + Vii-fD^ - M)tt, (2) 

where = d^Ay — d^A^, = <9 M — ieA^, A^(x) is an Abelian gauge field, ^ and ^ = ^+7° 
are two-component fermion fields of mass M, e is the coupling constant, and the matrices 7 M 
are chosen in the form 

In terms of Lorentz coordinates, the Lagrangian density for the boson theory equivalent to 
QED-2 can be written in the form [21] 

L=i(W-mV) + ^: e ^:+^:e^: ) 7 = m = (4) 

where C = 0.577216 is the Euler constant, 9 is a quantity that parametrizes the 9 vacuum of 
the fermion formulation of the theory, and the normal-ordering symbol means that diagrams 
with connected lines are excluded in perturbation theory in 7 (this corresponds to the usual 
meaning of the normal-ordering symbol in the Hamiltonian) - it is equivalent to perturbation 
theory in M. 

In [21, 22], the light-front Hamiltonian generating a theory that describes a one-component 
fermion field ip and which, in the limit where the regularization is removed, is equivalent in all 
orders in 7 to the Lorentz-covariant theory specified by the Lagrangian density (4) was found 
by using the method described in the Introduction. The theory defined by this Hamiltonian 
was regularized by the discrete-light-cone-quantization method mentioned in the Introduction: 
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the cutoff \x | < L and the antiperiodic boundary conditions in x were introduced for the 
field if>. The resulting corrected light-front Hamiltonian has the form 

H=J dx-(£ (dZ l [^}f - ^ ( e -W^) e-4 + h.c) - ^V^V), (5) 

—L 

where [. . .] denote the omission of the zero mode in x~ . The field if) is expanded in terms of 
creation and annihilation operators as 

iP(x) = -L f £ b n e-^-\)*- + Y: d^H)*- ) , (6) 

\n>l n>0 J 



{b n , &+} = R, d+} = <W, b n \0) = d n \0) = 0. (7) 
The operator u) is the quantity canonically conjugate to the charge operator Q 

Q = E K bn - E <dn, (8) 
n>l n>0 

which specifies the physical subspace of states, |phys): 

Qlphys) = 0. (9) 
The operator uj possesses the properties that completely define it [9, 19], 

e-|0) = &+|0>, e-- 10) = 4|0) (10) 

and 

e iw iP{x)e~ lw = e^Xx), (11) 

whence it follows that 

e^b n e-^ = b n+1 , e^d+e-™ = d+_ x , n > 1, e^d+e'^ = h. (12) 

The parameter 9 appearing in the Hamiltonian in (5) is a function of the ratio Mj e and the 
vacuum angle 9. This function is defined as a perturbation-theory series in M; therefore, its 
an explicit form remains unknown. Details concerning the appearance of the parameter 9 in 
the Hamiltonian are considered in the Appendix. Among other things, it is established there 
that, in the first-order in M, we have 9 = 9 and that, at any value of M, the parameter 9 is 
an odd function of 9 and takes the value of 9 = n at 9 = ir. In particular, it follows from the 
oddness of the function 9(9) that 9 = at 9 = 0. It should be noted that the parameter 9 can 
be related to the values of the vacuum condensates in the Lorentz-covariant theory [21, 22]. 

In calculating the mass spectrum of bound states, the quantity 9 is an independent pa- 
rameter of the theory, along with M and e; the relation between 9 and 9 for 9 ^ 0, ir can in 
principle be found by comparing the results obtained by calculating the mass spectrum of the 
theory in Lorentz coordinates and the theory on the light front. We note that expression (5) 
for the corrected light-front Hamiltonian and the expression obtained upon the naive canonical 
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quantization of the original fermion theory (2) in the light-front coordinates differ only by the 
addition of the second term, which is linear in the field operators and which depends on 6 and, 
hence, on the vacuum angle 6. Thus, we see that the naive canonical quantization does not 
take into account vacuum effects. 

The lightlike-momentum operator P_ has the form 

P = £ Kbn\ [n - i) + XX4| (n + \) . (13) 
This expression is used to calculate the mass spectrum of bound states. 



3. Calculation of the mass spectrum of bound states 

In order to find the mass spectrum of bound states of the theory, we will seek the eigenval- 
ues Ei of the fermion light-front Hamiltonian (5) (the subscript i numbers eigenvalues in the 
ascending order), 

H\*i) = Ei\*i) (14) 

in the subspace of physical states at a fixed value of the lightlike momentum (13). This subspace 
is specified by the conditions 

Q|tf)=0, P_|tf)=p_|tf). (15) 

In view of the antiperiodic boundary conditions and the first of the equalities in (15), the 
eigenvalue p_ has the form 

V- = ^N, (16) 
where N is a nonnegative integer. The bound-state masses Mj are given by 

Ml = 2p_Ei = Y NEi - ( 17 ) 

If, in expression (5) for the Hamiltonian, one performs the change of integration variable x~ = 
—z and uses expansion (6), the operator ^NH, which determines the quantities Mf does not 
involve the regularization parameter L explicitly, but it depends on N. The quantity L affects 
the mass spectrum only through relation (16). Since p_ does not depend on L, one can deduce 
from (16) that the limit iV — > oo corresponds to the limit where the regularization is removed, 
L — > oo. 

Since we use the antiperiodic boundary conditions, the subspace specified by the conditions 
in (15)) appears to be finite. This occurs because there exists a minimum positive value of the 
lightlike momentum P- — and because the creation operators correspond only to positive 
values of p_. An arbitrary state satisfying the conditions in (15) has the form 

\y) = d+_ 1 hN ...df N ^bf N ...b+ h \0), h = 0,l; 

2N / i N 2N j 



£sign iV-* + -U = 0, £ 
i=i v z/ i=i 



N-i +2 



k = N, (18) 
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where sgn is the sign function. The values U takes here only the values of or 1 by virtue of 
the anticommutation relations (7). 

The finiteness of the subspace that is specified by formulas (18) and where it is necessary to 
solve Eq. (14) reduces the problem to finding the eigenvalues of an iV ma t x iV ma t finite matrix 
whose elements are specified by the matrix elements of the operator ^NH between the states 
in (18). A precise solution to this problem can be found numerically. The resulting eigenvalues 
will determine the squares Mf of bound-state masses. It should be noted that the dimension of 
the matrix, JV mat , grows fast as the parameter N increases - this is reflected in Table 1. In this 



N 


9 


10 


12 


14 


16 


18 


20 


22 


24 


26 


28 


29 


30 


-A ma t 


31 


43 


78 


136 


232 


386 


628 


1003 


1576 


2437 


3719 


4566 


5605 



Table 1. Relation between the parameter N and the dimensionality N mSLt of the space of states. 

study, the maximum achieved values of N are N = 30 for the cases of 9 = 9 = and 9 = 9 = n 
and N = 28 for the remaining cases. 

As was mentioned above, the limit N — > oo corresponds to the removal of the regularization. 
For this reason, it is not sufficient to calculate the mass spectrum Mf at the maximum accessible 
value of N - it is necessary to analyze the behavior of the spectrum as a function of N and 
to find the way to extrapolate the calculated values to the region of N — > oo. We propose 
the following method of extrapolation. We introduce the quantity u — 1/N and consider 
the function M?(u). It is necessary to extrapolate the values of this function to zero. Our 
calculations reveal that the bound-state masses are sensitive to the parity of N - that is, there 
can occur a sharp change in Mf in response to the reversal of the parity of N. Therefore, it 
is reasonable to extrapolate the function Mf{u) to zero by two methods, individually in even 
and in odd values of N. 

In order to extrapolate the function Mf(u) to zero, we approximate the dimensionless ratio 
M 2 {u)/e 2 by polynomials of various degrees by the least squares method. We denote by Pj(ra) 
the value of a polynomial of degree n at zero. It is obvious that the maximum degree n that 
can be used is less by unity than the number of points at which the approximated function 
is known. In the present study, this degree is equal to 10 for the cases of 9 = 9 = and 
9 = 9 = ii and to 9 in the remaining cases (as the minimum value of N, we adopt iV = 9 in our 
calculations). 

At different values of the ratio M/e and the parameter 9, there arise different types of 
behavior of Pj(ra) as a function of n. In some cases, the function Pj(ra) tends to a saturation 
and changes slowly with increasing n (see Fig. la). In these cases, the value at which the 
saturation occurs will be considered as the result of extrapolating the function M?(u)/e 2 to 
zero. 

Sometimes, there arise oscillations against the background of a saturation (see Fig. 16). This 
occurs if the error with which one calculates Mf(u) becomes sizable. In this case, a function 
that involves a noticeable random noise is approximated by a polynomial of high degree. This 
usually takes place at large values of the ratio M/e, in which case the coefficients of various 
terms of the Hamiltonian (5) differ from one an other considerably. Conceptually, this situation 
does not differ from the preceding one. The value obtained by averaging the oscillations in the 
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n N 



Fig. 1. Examples of the dependence of the extrapolated value of P\(n) on the degree n of the 
approximating polynomial (a-c) and example of the dependence of the mass of the lowest bound 
state on N (of) for the following parameter values: (a) 9 = 9 = and M/e — 1, (b) 6 — — 
and M/e = 2 7 , and (c,d) 9 = 9 = it and M/e = 0.5. The symbols o and * represent the results 
of the extrapolation in, respectively, even and odd values of N. 



region of saturation will then be treated as the result of the extrapolation of the function 
Mj(u)/e 2 to zero. It is obvious that the value found in this way will have an error larger than 
that in the preceding situation. 

In the remaining cases, there arises the situation where the saturation cannot be seen - 
the function Pj(n) strongly changes with increasing n (see Fig. lc). In order to discriminate 
between this case and the two preceding versions of behavior of the function Pi(n), we introduce 
the measure of its relative variation in response to a considerable variation in the degree n of 
the polynomial (in the region of accessible values), for example, in the form 



^ ^(^(4)2 + ^(g) 2 ) /2' (19) 

This quantity characterizes the error with which the calculated values of the mass spectrum 
Mf describe its limiting value. We will assume that, if £ < 0.1, the saturation takes place for 
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the function Pi(n), so that its value obtained at accessible n describes well the limiting value 
of the bound-state mass. But if £ > 0.1, the calculated values of Pj(ra) do not characterize 
the behavior of Mf in the limit iV — > oo. Moreover, the most frequently occurring form of the 
dependence Pi(n) suggests its linearly decreasing character (Fig. lc). The same applies to the 
original dependence of the quantity Mf on the parameter N (see Fig. Id). This gives sufficient 
grounds to assume that, at these values of the ratio M/e and the parameter 9, the quantity 
Mf tends to — oo in the limit N — > oo. The possible reasons behind this effect are discussed in 
the next section. 

It should be noted that the choice of the value of 0.1 as a boundary one for the error in £ and 
the specific choice of formula (19) are arbitrary to a considerable extent, but, unfortunately, 
our calculations could not provide a more rigorous way to discriminate between the situations 
where the limit of the quantity Mf for N — > oo exists and where Mf tends to — oo in the same 
limit. 



4. Results of the calculations 

4.1. Case of 9 = 9 = 0. 

In the case of 9 = 0, the mass spectrum of the massive Schwinger model in Lorentz coordi- 
nates has received quite an adequate study (see [13, 16] and references therein). Usually, one 
studies the masses Mi and M 2 of the first two bound states, which are referred to as a vector 
and a scalar state, respectively. The most accurate results were obtained in [14] with the aid 
of lattice calculations. Table 2 presents the values of (Mi — 2M)/e and (M 2 — 2M)/e (it is 
precisely these quantities that were calculated in [14]) that were found by the method proposed 
here (with the aid of an extrapolation to the limit iV — > oo) and the values of these quantities 
from [14]. In the case of 9 = 0, the error £ does not exceed the threshold of 0.1 at any values 
of the ratio M/e that were considered here (more specifically, £ < 0.01 for Mi and £ < 0.03 for 
M 2 ) - that is, the above procedure of extrapolation to the limit iV — > oo provides quite reliable 
results. 

Figure 2a gives the extrapolated values of (Mi — 2M)/e along with the results obtained 
in [14]. The displayed errors were found from the corresponding values of the relative error 
£. In this figure, we also present the results corresponding to the maximum accessible value 
of iV = 30 - that is, the results obtained without extrapolation. It can be seen that these 
results are accurate only at small values of M/e; at the same time, the extrapolated values 
give a very good result up to M/e = 8. For M/e > 8, the extrapolated values reproduce a 
correct result within the error, which begins growing fast in this region. This is because the 
ratio Mj e becomes large in this region, with the result that the absolute error of the difference 
(Mi — 2M) / e appears to be large even at a small relative error in the calculated quantity Mf / e 2 . 
In order to depict the calculated results over the entire wide region of Mj e, it is convenient to 
plot the normalized values 

M norm = Mi (20) 

y/m? + (2M) 2 

as was proposed in [16]. These normalized values possess the property that M" orm — > 1 both 
in the limit M — > and in the limit e — > 0. In Fig. 2b, we show the curves for Mf orm that 
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Fig. 2. Calculated mass M 1 of the vector bound state at 9 = 9 = 0; * - results obtained 
by extrapolation to the limit iV — > oo, A - results of the calculation at N — 30, □ - results 
borrowed from [14]. 
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M/e 


(Mi - 2M) /e 


(M 2 - 2M) /e 


our study 


[141 

L J 


our study 


[141 

L J 


2 -io 


0.564 




1.13 






0.564 




1.13 




2-8 


0.563 




1.13 




2- 7 


0.563 




1.14 




2-6 


0.561 




1.15 




2- 5 


0.559 




1.17 




2-4 


0.554 




1.20 




2" 3 


0.545 


0.543 


1.23 


1.22 


2- 2 


0.524 


0.519 


1.24 


1.24 


2- 1 


0.489 


0.485 


1.20 


1.20 


2° 


0.445 


0.448 


1.12 


1.12 


2 1 


0.393 


0.394 


0.99 


1.00 


2 2 


0.339 


0.345 


0.84 


0.85 


2 3 


0.295 


0.295 


0.75 


0.68 


2 4 


0.279 


0.243 


0.74 


0.56 


2 & 


0.302 


0.198 


0.84 


0.45 


2 6 


0.368 




1.08 




2 7 


0.497 




1.41 




2 8 


0.619 









Table 2. Masses of the vector [M\) and scalar (M 2 ) bound states according to the calculations 
performed in the present study and in [14] at 9 = 9 = and various values of the ratio M/e. 



correspond to those in Fig. 2a. One can see that the extrapolated values agree with the results 
from [14] to a high precision over the entire range of Mj e. 

In Fig. 3, the extrapolated values of M" orm that were calculated on the basis of the corrected 
Hamiltonian (5) are contrasted against the analogous values corresponding to the Hamiltonian 
obtained upon the naive light-front canonical quantization of the fermion theory specified by 
Eq. (2) that is, to expression (5) where one discards the second term [see the comment before 
formula (13)] and those corresponding to the Hamiltonian obtained upon the naive light-front 
canonical quantization of the boson theory specified by Eq. (4) that is, to expression (5) where 
one discards the third term and replaces 9 by 9, as was shown in [21, 22]. One can see from this 
figure that the naive light-front canonical quantization of the fermion formulation of the theory 
gives good results at large values of the ratio Mj e - that is, in the region of weak coupling-while 
the analogous quantization of the boson formulation of the theories gives good results at small 
values of this ratio - that is, in the region of strong coupling. In order to obtain a light-front 
Hamiltonian that would provide good results at any values of the ratio M/e, it is necessary to 
implement the procedure of correcting the naive Hamiltonian, as was done in [21, 22]. 

In Fig. 4, we give the curves that represent the mass M 2 of the scalar bound state and which 
are analogous to those in Fig. 2a for the mass M 1 of the vector bound state. The behavior 
of the curves is identical to that in the case of the vector state: the values corresponding to 
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Fig. 3. Calculated mass M 1 of the vector bound state at 9 = 9 = 0; * - results obtained with 
the aid of the corrected light-front Hamiltonian, o and y - result obtained by means of the 
naive light-front quantization of, respectively, the fermion and the boson formulation of the 
theory. 




Fig. 4. Calculated mass M 2 of the scalar bound state at 9 = 9 = 0; * - results obtained by 
means of an extrapolation to the limit N — > oo, A - results corresponding to N = 30, □ - 
results borrowed from [14]. 
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N = 30 give good results only at small values of Mj e, while the extrapolated values give very 
good results up to M/e = 4 and reproduce the correct result within the error for M/e > 4. 
4.2. Case of 9 = 9 = vr. 

The value of 9 = it is of particular importance in the theory. It was predicted in [18] that, 
at 9 = re, a phase transition occurs in the theory at some value of the ratio M/e, so-called semi 
asymptotic fermions appearing in the region below the phase transition. In the region above 
the phase transition, as well as in the case where 9 ^ ir, confinement takes place. More recent 
calculations, performed for 9 = n revealed (see, for example, [13]) that the phase transition 
occurs at M/e = 0.33. 

By means of lattice calculations, the mass of the lowest state in the electron-positron (two- 
particle) sector as a function of M/e at 9 = tc was studied in [13]. Within our approach, this 
mass corresponds to the quantity Mi. 

In Fig. 5, the extrapolated values of M 1 /e are given along with the results reported in 
[13]. One can see that, at small values of the ratio M/e, these results agree well, but that 




Ql , , , , , , , L_ 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

M/e 



Fig. 5. Calculated mass M\ of the lowest bound state at 9 = 9 = ir; * - results obtained by 
means of an extrapolation to the limit iV — > oo, □ - results borrowed from [13]. 

this agreement deteriorates as the ratio M/e grows. For M/e > 0.14, the extrapolation error 
£ [see formula (19)] exceeds 0.1; therefore, the corresponding points on the graph are not 
quite reliable. For M/e > 0.21, there sharply appear very large oscillations in the behavior 
of the corresponding functions P\{n) describing the dependence of the extrapolated values 
on the degree of the approximating polynomial [see Section 3, the text before Eq. (19)]. At 
M/e = 0.4, the oscillations virtually disappear, whereupon the dependence begins displaying a 
linearly decreasing character. By way of example, this dependence at M/e = 0.5 is depicted in 
Fig. lc, while the corresponding dependence of the quantity Mf/e 2 on N (it is also manifestly 
linear) is shown in Fig. Id. For M/e > 0.2, the extrapolation error £ exceeds 0.5. 
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The linear decrease in M 2 /e 2 with increasing N gives sufficient grounds to assume that 
this quantity tends to — oo in the limit N — > oo. This means that, at the above values of the 
ratio M/e and the parameter 9, the spectrum of the Hamiltonian in (5) upon the removal of 
the regularization appears to be not bounded from below, so that the theory specified by this 
Hamiltonian is incorrect. This situation is possible in the case where there arise effects that are 
purely nonperturbative from the point of view of perturbation theory in the fermion mass M 
since the Hamiltonian in (5) was constructed by analyzing such a perturbation theory (in all 
orders). Obviously, the presence of the aforementioned phase transition is a nonperturbative 
effect in this case. One can conclude that above the phase-transition point (M/e = 0.33), the 
theory generated by the light-front Hamiltonian (5) becomes incorrect; at the same time, the 
original theory in Lorentz coordinates, which is specified by Eq. (2), remains correct, this being 
corroborated by the results reported in [13] for the region above the phase-transition point. 

The appearance of the aforementioned strong oscillations of the functions P\ (n) in the range 
0.2 < M/e < 0.4 is likely to be associated with the proximity of the phase-transition point, 
where the regularization, which is parametrized by the number N, can distort the theory more 
strongly than as usual. It can also be conjectured that a sizable deviation of the extrapolated 
values of M 1 /e from the results of paper [13] in the upper part of the region M/e < 0.2 is due 
to the same factor. 

4.3. Case of Intermediate Values of 9. 

As was indicated above, the quantity 9 is a function of the ratio M/e and the parameter 
9, this function being specified in the form of an infinite series in M. Therefore, the relation 
between 9 and 9 is a priori unknown [this is not so only in the particular cases of 9 = 0, n (see 
above)]. In principle, this relation can be sought by comparing the mass spectrum calculated on 
the basis of the light-front Hamiltonian (5), which depends on 9, and the spectrum calculated 
in Lorentz coordinates at a specific vacuum angle 9. 

It can be shown that the mass spectrum of the theory is invariant under the reversal of the 
sign of the quantities 9 and 9 (it should be recalled that 9 is an odd function 9). This can be 
seen most straightforwardly in the boson form of the theory (4), where the reversal of the sign 
of 9 is equivalent to the replacement (ip by — ip, which does not change the mass spectrum. 
However, this invariance can be directly seen from the Hamiltonian in (5). One can show that 
it does not change if we reverse the sign of 9 and simultaneously interchange the operators d n 
and b n+ i (in this case, u is replaced by —uo). This is a unitary transformation and does not 
change the mass spectrum. Thus, it is sufficient to perform calculations of the mass spectrum 
for the case where the parameter 9 lies between and n. 

For the lowest bound state, the mass Mi calculated in this way and normalized according 
to (20) is displayed in Fig. 6. Each curve corresponds to a fixed value of the quantity 9 from 
the set 0, 0.057T, 0.17T, . . . , n. Each successive curve (for increasing 9) lies below the preceding 
one. In the cases of 9/n = and 9/n = 0.05 (the first and the second curve from above), the 
extrapolation error £ [see (19)] does not exceed 0.1 at any values of M/e. 

In the case of 9/n — 0.1, £ > 0.1 at M/e values in the range between 2 and 2 4 ; therefore, 
it is meaningless to plot the corresponding points on the graph, so that the curves decompose 
into two parts. In this region, the corresponding functions P\{n) (used in the extrapolation) 
display a manifest linearly decreasing character. This situation is similar to that considered in 
Subsection 4.2 for the case of 9 = it. It can be concluded that, in the region being considered, 
the theory specified by the Hamiltonian in (5) becomes incorrect upon the removal of the 
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Fig. 6. Calculated mass Mi of the lowest bound state as a function of M/e at various values 
of the parameter 9: 0, 0.057T, 0.17T, . . . , ir. As the parameter 9 is increased, the curve is shifted 
downward. 



regularization, and it can be assumed that, in the vicinity of the point M/e = 2, there exists 
some nonperturbative effect as in the case of 9 = n (see Subsection 4.2). 

In principle, it can be assumed that a nonperturbative effect exists in the vicinity of the 
point M/e = 2 4 as well, above which the error £ again falls below 0.1. However, it seems more 
probable that, in the region of large values of M/e, the Hamiltonian remains unbounded from 
below upon the removal of the regularization, but the decrease in the mass of the lowest state 
with increasing regularization parameter N is so slow that £ appears to be less than 0.1. This 
is favored by the dependence Pi(n) which, at large values of M/e, has the form of a linear 
function with a moderate slope that yields £ < 0.1. As was discussed at the end of Section 3, 
there is unfortunately no method for discriminating between the situations where, in the limit 
iV — > oo, there exists the limit of the mass of the lowest state and where this mass tends to 
— oo. 

In the cases of 9/ir — 0.15, . . . , 0.3, the situation is perfectly analogous to that in the case 
of 9/n = 0.1 (only the width of the region where £ > 0.1 changes), which was considered 
immediately above, while, at 9/n = 0.35, . . . , 1 the difference consists in that the region where 
£ < 0.1 at large values of M/e is not reached at the M/e values considered here. 

In Fig. 7, the points where the extrapolation error £ is less than 0.1 are shown in the plane 
spanned by the parameters M/e and 9. It is natural to assume that, in the region where there 
are no points in the figure, the light-front Hamiltonian (5) is unbounded from below upon the 
removal of the regularization, N — > oo. 
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Fig. 7. Set of pairs of M/e and 6* values at which the extrapolation error £ is less than 0.1. 



5. Conclusion 

A numerical nonperturbative calculation of the mass spectrum of QED-2 (massive Schwinger 
model) has been performed by the method of discrete light-cone quantization. In doing this, 
use has been made of the corrected light-front Hamiltonian (that is, that which generates a 
theory that is equivalent in all orders of perturbation theory in the fermion mass M to the 
Lorentz-covariant formulation of QED-2) constructed in [21, 22]. The results obtained in this 
way have been compared with the results of the numerical calculations on a lattice in Lorentz 
coordinates from [14, 13]. 

Since actual calculations are performed at a finite value of the infrared-regularization pa- 
rameter N, a method has been proposed for extrapolating the results of these calculations to 
the region of N tending to infinity, this corresponding to the removal of this regularization. 
As the result of this extrapolation, it becomes clear that, at some values of the parameters of 
the theory, its mass spectrum is not bounded from below in the limit where the regularization 
is removed. This occurs only at rather large values of the fermion mass M; therefore, it is 
natural to assume that, in this region, the Hamiltonian used becomes incorrect, since it was 
constructed via an analysis of perturbation theory in M. 

The calculations have been performed over a wide range of the fermion mass M for all values 
of the Hamiltonian parameter 9, which is a function of the ratio M/e and the vacuum angle 0. 
This function, which is a priori unknown possesses the property that it vanishes at = and 
is equal to n at = it. 

At zero value of the vacuum angle 0, the resulting spectrum is bounded from below at any 
value of M, and the values obtained in this case for the two lowest bound states reproduce well 
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the results reported in [14]. 

At the vacuum-angle value of 9 = n, which is special for the theory being considered, the 
masses found here for the lowest bound state agree well with the results of paper [13] for rather 
low values of M; as the mass M increases, there first arises a discrepancy, whereupon the 
spectrum of the theory become unbounded from below. Since the M value at which this takes 
place is approximately equal to that at which there occurs a phase transition in the theory 
(see, for example, [13])), it would be reasonable to assume that, at M values above the phase- 
transition point, the light-front Hamiltonian used here, which was constructed on the basis of 
an analysis of perturbation theory in M, becomes incorrect. 

Our calculations have revealed that, in the case of QED-2, the procedure employed in 
[21, 22] to construct the corrected light-front Hamiltonian leads to a Hamiltonian that yields 
good results in nonperturbative calculations. By exploring the question of how the spectrum 
of the Hamiltonian changes in response to changes in the parameters of the theory from the 
perturbative region of their values, one can determine the boundaries of the applicability region 
of this Hamiltonian - that is, to find the region where it is necessary to take additionally into 
account nonperturbative (for example, vacuum) effects. This may be of use in studying more 
realistic gauge field theories. 

The LF Hamiltonian (5), used in the present paper for the calculation of bound state mass, 
includes the operator u, which has no simple expression in terms of field operators. It is defined 
only by its properties (10)-(12). Due to this fact the expression for the Hamiltonian depends 
essentially on the form of the regularization, i.e. \x~\ < L and antiperiodic boundary conditions 
in x~ for the field ip. Now we have found a possibility to rewrite the expression (5) in such a 
way that it contains only fermion field operators, and describes in the limit of removing the 
regularizations the same theory as the Hamiltonian (5). This new expression has at 9 = 9 = 
the following form: 

H=f dx- (dZ 1 ^}) 2 + ^ (d+do + bfh) - ^V^lV) • (21) 

Preliminary calculations of the mass spectrum, produced by this Hamiltonian, show that in the 
limit of removing the regularization, L — > oo, results indeed coinside, with a good accuracy, 
with those for the bound state mass spectrum found here for the Hamiltonian (5). Work on 
studing of the Hamiltonian (21) spectrum and also on the constructing the analogous expression 
for the Hamiltonian at 9 ^ will be continued in future. 
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the Russian Foundation for Basic Research, grant no. 05-02-17477 (S.A.P. and E.V.P.). 

Appendix 

Let us find out how the parameter 9, which is a function of the ratio M/e and the vacuum 
angle 9, appears in the light-front Hamiltonian (5). It was shown in [21, 22] that the coefficient 
of the operator e^d^ in the integrand on the right-hand side of (5) is the limit of the quantity 
— B* (* denotes the operation of complex conjugation) in the limit w — > oo (which corresponds 
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to removing the intermediate ultraviolet regularization), where B is given by 




1 + ^_A" 2 + iA". (A.l) 



Aw 2 w 

Here, A' and A" are, respectively, the real and the imaginary part of the sum (calculated in 
Lorentz coordinates and in all orders in 7, including the first order) 

A = 7^ + £^7 fe (A.2) 

1 k=2 

of all connected diagrams of the boson theory specified by the Lagrangian (4) which have a 
property that all their external lines connect with the one vertex ^e ld [the propagators of the 
external lines and a common factor that depends on their number are not included in the 
definition of A, as can be seen from the first term in (A.2)]. In [21], it was established that, in 
the limit w — > 00, the quantities A k for k 7^ 2 are finite, while A 2 behaves as 

7 2 

A 2 = —w + const. (A. 3) 

4 

Expression (A.l) was deduced from an analysis of perturbation theory in 7 (in all orders). 
As a matter of fact, the series that can be obtained by substituting expansion (A.2) into (A.l) 
provides the definition of the quantity B. One can see that the first term of this series is linear 
in 7 and, at large w, its radius of convergence varies in proportion to \ jw. As was mentioned 
above, it is necessary to find the limit of the quantity B for w — > 00. Since the aforementioned 
radius of convergence tends to zero in this limit, it is obvious that, before going to the limit, 
the quantity B must be continued analytically in 7 to the region of positive 7 values, which lie 
beyond the disk determined by this radius of convergence. For this, we determine the behavior 
of radicand on the right-hand side of (A.l) at large w and 7 of about l/w. Taking into account 
(A. 3), we obtain 

/ 1 A I 1 / / 1 \ 

(A.4) 

From here, we find that there exist two branch points, whose positions are given by the formula 

71,2 = —n+o(^H). (A.5) 

wcos9 \w 6 l 2 J 

It can be concluded from formulas (A.4) and (A.5) that, in the limit w — > 00, the sought 
analytic continuation of the quantity B has the form 




B = sign(cos A" 2 + iA" = ^e i§ . (A.6) 

The form of this expression corresponds to the form of the coefficient of the operator e^dg 
in Hamiltonian (5). Expression (A.6) differs by the presence of the sign function sign(cos#) 
from the corresponding expression presented in [21, 22], where the features of the analytic 
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continuation of B were not taken into account. In the first order in 7 (and, hence, in M), we 
find from (A. 6), with the aid of expansion (A. 2), that 9 = 9. 

Upon the Euclidean rotation in the diagrams determining the quantity A, it becomes clear 
A is a real- valued function of m, 7, e %e and e~ %e . It follows from here that the reversal of the sign 
of the parameter 9 is equivalent to the complex conjugation of the quantity A, this quantity 
being real- valued at 9 = ir. By using these facts and formula (A. 6), one can easily see that 9 is 
an odd function of 9 and, in addition, that 9 = n at 9 = it. 
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